The energy dependence of the pp — > K + nT, + reaction close to threshold 
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The production of the E + hyperon through the pp — > K + nT, + reaction has been investigated 
at four energies close to threshold, 1.826, 1.920, 1.958, and 2.020 GeV. At low energies, correlated 
K + tt + pairs can only originate from E + production so that their measurement allows the total 
cross section for the reaction to be determined. The results obtained are completely consistent 
with the values extracted from the study of the _K" + -proton correlation spectra obtained in the 
same experiment. These spectra, as well as the inclusive K + momentum distributions, also provide 
conservative upper limits on the E + production rates. The measurements show a E + production 
cross section that varies roughly like phase space and, in particular, none of the three experimental 
approaches used supports the anomalously high near-threshold pp — > K + nT> + total cross section 
previously reported [T. Rozek et al, Phys. Lett. B 643, 251 (2006)]. 



PACS numbers: 13.75.-n, 14.20.Jn, 14.40. Aq, 25.40.Ve 

I. INTRODUCTION 

The energy dependence of the total cross section for 
associated strangeness production in the exclusive pp —> 
K + pA reaction near threshold has been investigated in a 
series of experiments undertaken at the COSY accelera- 
tor of the Forschungszentrum Jiilich using the COSY-11 
and COSY-TOF detectors P-Q. The resulting behavior 
can be described in terms of a three-body phase space 
that is influenced by a strong final state interaction (FSI) 
between the A hyperon and the emerging proton. 

The situation for E° production in the analogous pp — > 
K + pT,° reaction is rather different since the existing data 
show a total cross section that varies as expected from 
three-body phase space, with no evidence for a strong 
T,°p FSI [H, [H, 0] • A second distinction between the pro- 
duction of the two hyperons in proton-proton collisions 
is that, at the same value of the excess energy e close 
to threshold, the total cross section for pp — > K + pYP 
is at least one order of magnitude less than that of 
pp — » K + pA. This probably reflects differences in the 
underlying reaction mechanisms as well as between the 
two FSI involved. 
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Since the S has isospin 7 = 1, extra light may be 
cast on both the FSI and reaction mechanism questions 
through the study of S + production in the pp — > K + nT, + 
reaction. The first measurements on this channel were 
completed by the COSY-11 collaboration at proton beam 
energies of T p = 1.826 GeV and 1.958 GeV, i.e., at 
e = 13 MeV and 60 MeV, by detecting the neutron in 
coincidence with the produced K + Q. The values ob- 
tained for the total cross section were strikingly high com- 
pared to those for E° production, or indeed with those 
for A production, with ratios of i?(E + /S°) = <r(pp — > 
K+nJ: + )/a(pp -> K+pYP) = 230±70 at e = 13 MeV 
and 90±40 at 60 MeV @. Such values would imply a 
very anomalous energy dependence, since bubble cham- 
ber data yield a ratio of R(E + /S°) = 2.3 ± 0.9, though 
at the much higher excess energy of e = 350 MeV @. 

In an effort to understand the situation, a first mea- 
surement of the pp — > K + nT, + reaction was carried out 
at the COSY-ANKE facility that did not rely on the 
detection of the final neutron 0j. This experiment at 
e = 129 MeV (T p = 2.16 GeV) involved a three pronged 
approach. Although the statistics were low, apart from 
a small K + nAir + background, the detected K + ir + co- 
incidences at this energy could only have come from 
K + nTi + , where the hyperon decays through S + — » nir + 
(BR = 48.3%). This approach indicated that i?(E+/S°) 
was of the order of unity at this energy. 

Information on the K + nS + channel could also be ob- 
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tained through the measurement of AT + -proton corre- 
lations, where the E + shows up as a contribution in 
the missing-mass distribution via its decay E + — > pir° 
(BR = 51.6%). Although the overall spectra are sensitive 
to the model used to describe the dominant K+pA chan- 
nel, at high K+p missing masses there can only be contri- 
butions from £ production. Both this study and that of 
the inclusive K + production showed consistency with the 
K +n+ determination and clearly excluded a large value 
of i?(E+/E°). 

These first ANKE results were obtained at a signifi- 
cantly higher excess energy than those of COSY-11 Q 
and so did not rule out the possibility of a strong thresh- 
old anomaly. The aim of the present work is to extend 
the earlier ANKE methodology by making four measure- 
ments that cover the near-threshold region as well as one 
below the E threshold to obtain data where the K + pA 
channel is the only source of associated strangeness pro- 
duction [lol ]. The beam energy range covered was from 
1.775 to 2.020 GeV. 

We start by summarizing in Sec. [II] the existing near- 
threshold data on hyperon production in proton-proton 
collisions and their interpretation in terms of phenomeno- 
logical models. The ANKE experimental facility used 
for the study of the pp — > K+n'E+ reaction and the tech- 
niques involved there are described in Sec. lIIIl Section lTVl 
presents first the results of the measurements of inclusive 
K+ production. The comparison of these data taken just 
above with those from just below the £ thresholds allows 
us to set an upper limit on the £ + production cross sec- 
tion at e — 13 MeV that is almost two orders of mag- 
nitude smaller than the COSY-11 result 0. The K+p 
missing-mass spectra lead to values of the total cross sec- 
tions for both A and £° production that are in agreement 
with the world data. The high ends of these spectra can 
only be populated by events arising from E production 
and these allow first estimates to be made of the E + cross 
section as well as conservative upper limits. Although the 
statistics here are reasonable, it is hard to quantify the 
systematic uncertainties and, in this respect, the most re- 
liable determination of the E + production cross section 
is through the study of K+ir+ coincidences which, at the 
energies of this experiment, can only arise from E + decay. 
The values obtained from this study are completely con- 
sistent with those found from the K+p coincidence data 
and fall below the upper bounds set by the inclusive data. 

The results for the total cross sections are discussed in 
Sec.fVl where they are compared to other measurements 
of A and E° production as well as those of E + . The 
clean determinations of the pp — > K+nYj+ cross section 
from the K+tt+ and K+p measurements give values that 
are more than two orders of magnitude lower than those 
reported in the literature through the detection of the 
neutron Q. Our conclusions and outlook for the future 
are presented in Sec. EH 



II. HYPERON PRODUCTION IN 
PROTON-PROTON COLLISIONS 

In the low energy regime (around T p s» 2 GeV) there 
are three associated-strangeness production channels al- 
lowed in proton-proton collisions, viz. K+p A, K+pT, , 
and K + nYj + . Measurements of inclusive K+ spectra 
lead to distributions where there are contributions from 
all three final states. There is no possibility of distin- 
guishing between the different reaction channels without 
modeling the distributions from all of them |ll|, |l2[ , and 
this can lead to severe ambiguities. The basic problem 
here is that the cross section for A production well away 
from threshold is much larger than those for the E close 
to threshold and there is no reason to believe that the 
associated distributions should vary like phase space. 

The conditions improve when two particles are de- 
tected in the final state. Signals from the pp — > K+pA 
and pp — > K+pYP reactions are usually identified through 
the corresponding peaks in the K+p missing-mass spec- 
tra 0] . The physical backgrounds in such measurements 
are mostly connected with the decays of A and E° hy- 
perons, which lead to a second proton in the final state 
(see Table [TJ) . The remaining part of K+p missing-mass 
spectra can be associated with the pp — > K+nYj+ reac- 
tion channel, where the E + — > pir° decay gives rise to the 
proton that is detected. It is important to note that the 
missing-mass range that is obtained using a proton from 
A decay does not extend as far as that linked to E pro- 
duction so that the study of the maximum missing-mass 
region allows one to derive values for the weighted sum of 
E° and E + production. Nevertheless, the contributions 
are rather small and systematic effects can be significant 
near this kinematic limit. 

However, below the threshold for pp — > K+nAir° pro- 
duction (T p — 1.975 GeV), the only source of the K+tt+ 
correlations is E + production and this provides a much 
more reliable way to identify the channel than the alter- 
native of detecting the neutron Q . Although the mea- 
surement of K+n+ coincidences provides a very clean sig- 
nal for the identification of the pp — > K+nT,+ reaction, 
the statistics achievable at ANKE are quite low 0, [HI ■ 



Final state 


K+pA 


K+pY? 


K+nT,+ 


T thr [GeV] 
BR(K+p) {%] 
BR{K+tv+) [%} 


1.582 
63.9 


1.794 
63.9 


1.789 
51.6 
48.3 



TABLE I: Characteristics of the three hyperon production 
channels in proton-proton collisions studied in this experi- 
ment. In addition to the threshold beam energies, also shown 
are the fractions of the cross sections (BR) associated with 
the K+-K+ and indirect K+p coincidences. 

It is evident from the inclusive K+ production data 
that there is a stron g Ap final state interaction at low 
invariant masses [ill; E3- This is reflected in the en- 
ergy dependence of the pp — > K+pA total cross section 
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analyzed in Refs. [Kg, |l4j, where the near-threshold ex- 
perimental data are enhanced compared to phase space. 
Significant influence of the Ap final state interaction has 
also been found in the analysis of pp — > K + pA Dalitz 
plots 0, |||. The measured distributions have been de- 
scribed using a A production model which, in addition 
to the strong FSI, includes the iV*(1650)-isobar in the 
production mechanism. It is believed that this is the 
dominant resonance in the ~ 2 GeV region 

The ANKE acceptance for the pp — > K + pA reaction 
changes significantly over the energy range where the 
data were collected. In addition, different models have 
been used in this analysis for the estimation of the A 
total cross section at different energies. A three-body 
phase space, modified by the Ap FSI, has been used 
for T p < 1.9 GeV. Above 1.9 GeV, a model for the 
pp — > K + pA reaction, similar to the one used for the 
2.16 GeV data 0], has been developed. This assumes the 
dominance of the 7V*(1650)-resonance in the production 
mechanism and a strong Ap final state interaction. In ad- 
dition to these effects, there is some anisotropy in the an- 
gular spectra of the kaon and direct proton in the K + pA 
final state (ll^- fljl . Unfortunately, unlike the case for the 
2.16 GeV data [18|. there are no measurements of the 
angular distributions directly at the energies of interest. 
The relevant parameters have therefore been determined 
from linear interpolation of the results of Refs. [l5l - fl8j 
at nearby energies. More details of the model and the 
parameters used in the Monte Carlo simulations can be 
found in Ref. [ijj. 

Extensive measurements of the total cross section for 
E° production close to threshold were carried out by the 
COSY-11 collaboration d, Some studies of differ- 
ential observables were undertaken at the COSY-TOF 
facility [2| , though most of the results have only been re- 
ported in theses [13, HE]. There is little evidence from 
these data that in the vicinity of threshold there are sig- 
nificant deviations from phase space. We therefore use a 
phase-space description when evaluating the acceptance 
for E° production at ANKE. Since there is no experi- 
mental information at all on differential quantities for 
the pp — > K + nH + reaction channel at low energies, a 
similar three-body phase-space description has been em- 
ployed for E + production. 

The pp — > AT°pE + reaction was studied at s = 126 [l9| 
and 161 MeV [20] during the COSY-TOF pentaquark 
searches. The value of the total cross section deduced 
from these data at the lower energy is quite close to that 
found for the K + pH° final state [l8[ and this already 
allows simple bounds to be placed on the cross section 
for the pp — > K + nS + reaction. 

There are only two independent isospin amplitudes for 
E production in proton-proton collisions and these can 
be taken to correspond to the I = 1/2 and 7 = 3/2 
combinations of the final EV pair. There is therefore a 
linear relation between the amplitudes for the production 



of the three possible final states: 

f(pp ->■ AT+nE+) + f{pp ->■ K Q pY 1 +) 
+ V2f(pp^K+pY,°) = 0. 



(1) 



This leads to a triangle inequality between the total cross 
sections @ , which can provide model-independent limits 
on the ratio of E + to E° production: 



y/aipp^ AT°pE + ) - y/2(r(pp K+pYP) 
<a(pp^ K + nT, + ) 



< 



yJo{pp ->• AT°pE+) + V 2 <r(pp -> K+pT, ) 



(2) 



The COSY-TOF result at e = 126 MeV [T^ suggests 
that at the energy of the earlier ANKE experiment at 
2.16 GeV the ratio i?(E+/E°) of the total cross sec- 
tions for E + to E° production should lie in the range 
1/6 < i?(E+/E°) < 6. Though weak, this bound is 
an order of magnitude smaller than that found by the 
COSY-11 collaboration, though at lower energies @. A 
direct measurement of the pp — > AT + rtE + cross section 
at this energy using the COSY-TOF spectrometer 21 1 
satisfies well this inequality, whereas an earlier one at 
e = 94 MeV seems to be much too large [22| . 



III. EXPERIMENT 

The experiment was carried out at the ANKE facil- 
ity [23l ] with an unpolarized proton beam provided by 
the Cooler Synchrotron and storage ring COSY (24[. Of 
the five beam energies used, T p = 1.775, 1.826, 1.920, 
1.958 and 2.020 GeV, the first, just below the threshold 
for the production of the E hyperon, was used for back- 
ground studies. The values of the excess energies e for 
the three final states are given in Table [TTJ 



Final state 


K + pA\K + pE u 


K + nY, + 


P 


T P 




e 




GeV/c 


GeV 




MeV 




2.546 


1.775 


70 






2.600 


1.826 


88 


11 


13 


2.700 


1.920 


122 


45 


47 


2.740 


1.958 


135 


58 


60 


2.806 


2.020 


157 


80 


82 



TABLE II: Values of the excess energies e for the three re- 
action channels at the five beam momenta p and energies T p 
used in this experiment. 

The ANKE magnetic spectrometer, located at an inter- 
nal target position of COSY, is composed of three dipole 
magnets. Dl and D3 bend the direct proton beam from 
the undisturbed COSY orbit to the ANKE target and 
return it back, respectively. The analyzing magnet D2 7 
placed between D 1 and -D3, was operated in this exper- 
iment with a field strength of B = 1.57 T. This magnet 
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FIG. 1: (Color online) Sketch of the relevant parts of the ANKE detector system showing the positions of the two bending 
magnets Dl and Z)3 and the target before the analyzing magnet D2. Only information from the forward (Fd) and positive 
side (Pd) detectors, including the stop counters in the side wall, was used in this experiment. 



deflects positively charged particles with laboratory po- 
lar angles up to 20° towards the positive side (Pd) and 
forward (Fd) detector systems. The high density of the 
unpolarized H 2 cluster-jet target [25| gave a typical tar- 
get thickness of ~ 10 14 atoms/cm 2 and this resulted in 
an average luminosity of the order of ~ 2 nb _1 s _1 . In 
the present analysis only the positive side (for the 7r + , 
K + and p identification) and forward detector systems 
(for proton identification) have been used. A sketch of 
the ANKE detector system is presented in Fig. m 

The positive side detector system consists of 23 start 
(Sa) and 21 stop (So) counters for time-of- flight (TOF) 
measurements and two multiwirc proportional chambers 
(MWPC) for track reconstruction and background sup- 
pression. Depending on the type and momentum of the 
particle, the tracking efficiencies of the Pd MWPCs are 
typically ~ 90%, the values being determined using ex- 
perimental data. The momentum resolution in the posi- 
tive side detector system is about 2% (FWHM). 

The first fifteen stop counters are part of the range tele- 
scope system designed for the identification of low mo- 
mentum K+ (p£ ax « 620 MeV/c for B = 1.57 T) [Hj]. 
The remaining six stop counters, joined together to form 
the side- wall detector (p max pa 920 McV/c for B 
1 .57 T), were used for n + and p identification on the 
basis of time-of-flight criteria. 

Each of the range telescopes is made up of a stop 
counter, an energy loss counter (AE), a so-called delayed 
veto counter, and two passive degraders made of copper. 
The thickness of the first degrader is chosen such that a 
K + deposits all its energy in the AE counter and stops 
either at the edge of it or in the second degrader. The 
products from the K + decay are registered in the delayed 
veto counter with the characteristic decay time of 12.4 ns. 
All the range telescopes, placed in the focal plane of the 
D2 magnet, are optimized for the measurement of a given 
K + momentum. The delayed veto criteria leads to a sup- 
pression of better than 10 -5 in the non-kaon background 
for both inclusive and coincidence measurements. The 



details of the K + identification using the delayed veto 
technique can be found in Ref. [2(| . 

Due to the unchanged D2 magnetic field, the geome- 
try of the positive side detector system was fixed with 
respect to D2 so that the acceptance of each telescope 
remained the same at all the beam energies where exper- 
imental data were collected. As a consequence, the same 
set of detector efficiencies and acceptances could be used 
in the analysis of the inclusive data collected at different 
energies. 

The ANKE forward detector system used for the K + p 
correlation measurements and the luminosity determi- 
nation consists of two multiwire proportional chambers, 
one drift chamber, and a hodoscope of scintillator coun- 
ters (27J. The tracking efficiency of the Fd chambers, as 



determined from the experimental data, is homogeneous 
over the Fd acceptance and is better than 95% for pro- 
tons. The ensemble of chambers used in the experiment 
allows the reconstruction of the proton momentum with 
a precision of ~ 2% (FWHM). The Fd hodoscope com- 
prises two layers of plastic scintillator counters shifted 
with respect to each other by half the width of a counter. 
The two-layer structure allows the time from individual 
counters to be calibrated with respect to each other by 
using tracks of particles which cross counters that over- 
lap in different layers. The time signals from individual 
forward counters are therefore aligned such that in the 
analysis the full hodoscope can be treated as a single 
counter [27) . 

In order to normalize the experimental data, and hence 
extract absolute cross sections for S + production, the 
proton-proton elastic scattering rate was measured us- 
ing the Fd detector. Due to the high count rate from 
this reaction channel, a dedicated trigger, prescaled by a 
factor of a thousand, was used. Having determined the 
momentum of a proton candidate in the Fd, the pp — > pp 
reaction was identified via the missing-mass technique. 
The elastic scattering peak, with a width of 50 MeV/c 2 
(FWHM) , is well separated from the inelastic background 
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arising from pion production. The numbers of protons 
scattered between laboratory angles of 6° and 9° were 
determined from a fit to the data by a Gaussian with 
polynomial background. Predictions for the pp elastic 
differential cross section from the SAID analysis |35j were 
then used to evaluate the luminosity. For all the beam 
energies in the experiment, the shape of the differential 
cross section as a function of laboratory angle agreed with 
the predictions to better than 3%. The luminosity was 
determined with an overall uncertainty of 7%. This takes 
into account a 5% estimate of the uncertainty from the 
SAID program, 3% uncertainty in the analysis algorithm, 
and 3% uncertainty in the Fd MWPC efficiency determi- 
nation. The values of the total luminosities accumulated 
are reported in Table IIVI 

The Runge-Kutta method of momentum reconstruc- 
tion was employed for the final analysis. This uses all 
the information from the chambers and the measured 
D2 magnet field map [28[ . The precise knowledge of the 
particle trajectories provided by this method allows one 
to calibrate absolutely the times between different start- 
stop counters in the Pd using experimental data with a 
7r + or any other cleanly selected particle. The internal 
delays for time signals do not change during the experi- 
ment and the required constants could be calculated us- 
ing time information from individual counters and corre- 
sponding tracks. After performing the time calibration, 
the invariant mass of a particle with a given track could 
be calculated. 

The time-of- flight spectra measured at 1.920 GeV for 
K + selected using the delayed veto technique are pre- 
sented in Fig. [2] summed over all stop counters. 

For the extraction of the inclusive double-differential 
cross section, the residual background was subtracted 
from the time-of-flight spectra measured in individual 
stop counters. The number of K + detected in an in- 
dividual telescope, N K + , has been determined using a fit 
of the TOF spectra that assumes a Gaussian peak from 
the kaons and a flat background. The background con- 
tribution to the individual stop counter distributions is 
less than 5%. 

The laboratory cross section has been evaluated for 
every momentum bin from 



dfl dp 



(T P ) 



N 



K+ 



Ap Afi L t ot e K + ' 



(3) 



where Af2 is the solid angle integration region and L to t 
the integrated luminosity. A momentum bin of Ap = 
24 MeV/c was imposed in the course of the data analysis. 
The efficiency of K + identification, t K + , is estimated on 
the basis of 



_tel scint MWPC w ^.acc 
c X 6 X fc X c 



(4) 



Due mainly to geometric acceptance factors, the 
delayed-veto technique leads to a telescope efficiency e tcl 
of about 30% for low momentum kaons but only 10% for 
the highest momenta (26|. The value for each telescope 




20 
TOF [ch] 



FIG. 2: Time-of-flight (TOF) spectra, measured at 1.920 GeV 
with a delayed veto, presented on a logarithmic scale. The 
data, measured in units of 44 ps per channel, have been cen- 
tered around zero. The hatched histogram shows the distri- 
butions of K + events for which there were K + p coincidences. 



has been determined with a precision of 5% by using K + p 
correlations recorded with a dedicated trigger. The scin- 
tillator counter efficiency e scmt is higher than 98% while 
that of the MWPC, e MWPC , varies between 90-95%, de- 
pending on the K + momentum. These efficiencies have 
been determined using experimental data. The accep- 
tance correction factor e acc has been estimated using a 
GEANT4 [H model of the ANKE spectrometer. The 
overall K + detection efficiency e K + in this measurement 
is independent of beam energy. 

The K + p missing-mass spectra, which depend on the 
data measured with the Fd detector, have been used to 
extract cross sections for the A and S° production chan- 
nels. For the selection of K + p pairs detected in the for- 
ward and positive detectors, the time differences between 
the K + in the Pd and the proton in the Fd have been 
calculated. A 3a cut has been applied to subtract the 
residual background from accidental coincidences. The 
background in the K + p coincidence spectra measured 
with the delayed veto is less than 2%. 

The K + time-of-flight spectrum measured in coinci- 
dence with the proton detected in the Fd at 1.920 GeV 
is shown in Fig. [5] by the hatched histogram. The overall 
reduction between the K + TOF spectra measured with 
and without the proton being detected in the Fd is due 
to the acceptance of this detector. The numbers of A and 
S° events, N ev , in the spectra were estimated using fits 
of the simulations to the K + p missing-mass spectra. The 
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FIG. 3: (Color online) Identification of two-particle events in the positive side detector at 1.920 GeV, where a cut of 9 K + < 12° 
has been imposed on the K + angle, (a) The difference At between measured and calculated times of the correlated particles. 
The selection of double coincidences in Pd is shown by the hatched histogram for K + p and the filled histogram (red) for K + ix + 
pairs, (b) The K + p missing-mass spectra measured in Pd. The hatched histogram represents K + p missing mass selected using 
the cuts shown in panel (a), (c) The invariant mass of a detected particle before and after (filled histogram) the selection of 
K + tv + coincidences using the criteria in panel (a). 



total cross sections were then evaluated from 



f total 



-ktotai e s 



(5) 



where L to t is the integrated luminosity, g sum the over- 
all efficiency that includes the detector efficiencies, the 
in-flight K + decay correction, and the total acceptance 
within a particular reaction model, as discussed in Sec.llll 

The rate of K + ir + coincidences is very small. There- 
fore, the K + p correlations detected in the positive detec- 
tor have been used to optimize the cut parameters and 
selection criteria. The K + and p tracks were selected 
using time-of-fiight criteria. In Fig. [3^ the difference be- 
tween calculated and measured time in the So counters 
is shown, with its center shifted to zero. Also presented 
in panels b and c are the distributions in missing mass, 
assuming that tracks are from K + p coincidences, and the 
invariant mass of the particles whose tracks are detected 
in the Pd. The hatched histogram in Fig. \Sjp represents 
K + p coincidences selected using a 3er cut on the time 
difference shown in panel a. The number of A and E° 
events, estimated using the clear peaks in the missing- 
mass distribution, hardly change. This indicates that 
the efficiency of the cuts is better then ~ 97 %. 

When the same cuts on the time difference are applied 
to the K + tt + candidates, the results found are shown 
in Figs. [3Jd and c by filled histograms. Since it is im- 
possible to get correct time information for events where 
both particles are detected in the same counter, these are 
dropped from the analysis. This effect is most serious for 
the K + tt + coincidences at lowest energy but even there 
it does not exceed 25%. An additional correction factor 
has been introduced to compensate for the loss. 

The number of K + ir + correlated pairs has been used 
to evaluate the total cross section for E + production on 



the basis of Eq. (|3J|. This includes a correction for the 
corresponding branching ratio (48.3%) reported in Ta- 
ble |U The acceptance calculation, as well as the reaction 
models used, are discussed in Sec. [XT] and given in more 
detail in Ref. [HI. 



IV. EXPERIMENTAL RESULTS 
A. Inclusive counting rates 

A first limit on the cross section for E + production 
can be obtained by comparing the inclusive K + count- 
ing rates just above the E thresholds with those just be- 
low. Figure [4] shows the ratio of the data obtained at 
1.826 GeV to those at 1.775 GeV as a function of the kaon 
laboratory momentum. In this way any questions regard- 
ing telescope efficiencies and acceptances are avoided. At 
1.826 GeV the K + mesons associated with E + produc- 
tion have a maximum laboratory angle of ??^+ x = 12.2° 
and a momentum range from 352 to 643 MeV/c. Essen- 
tially the whole phase space for E production is therefore 
sampled in the data, though the vertical angular limita- 
tion and, to some extent, the horizontal angle cut of 12° 
and the ±12 MeV/c momentum cut for each telescope re- 
duces the acceptance to about 23%. The simulations dis- 
cussed in Sec. |TT] suggest that the contribution to the ratio 
from A production should be essentially constant over the 
whole momentum range, independent of any assumptions 
made regarding the kaon angular distributions. 

The ratio of differential cross sections can be expressed 
by Eq. ©: 



R 



1.826 
1.775 



Al.826 . Ej.826 + E+ 
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Ai 
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Ax. 



775 
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The shape expected for the contribution from the com- 
bined pp — > K + pTP and pp — > K + nY* + channels is in- 
dicated, using total cross sections cr(S°) = = 
0.021 pb, cri.826(A) = 7.9 ^ and cti. 775 (A) = 6.1 pb. De- 
pending upon how the A level is drawn, this comparison 
gives a cross section ratio of (a(pp — > K + pYf > ) + u(pp — > 
K+nY>+))/<j{pp -> K+pA) = (4 ±2) x 10~ 3 , from which 
we deduce that the total cross section for E + produc- 
tion at 1.826 GeV (e = 13 MeV) is below 45 nb at the 
98% confidence level. This upper bound already excludes 
the COSY- 11 result Q by a very large margin. Even if 
the ratio of the S + /E° production cross sections were 
a factor of six, as allowed by the triangle inequality at 
e = 129 MeV, the resulting distribution would greatly 
overestimate the data shown in Fig. 01 
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FIG. 4: (Color online) Ratio of normalized count rates for 
inclusive K + production at 1.826 GeV to those at 1.775 GeV 
as a function of the momenta in the different telescopes. The 
data were taken using the ANKE range telescope system with 
■& K + < 12° and a 12 MeV/c momentum cut for each tele- 
scope. The black solid histogram represents the simulation 
for the pp — > K + pA reaction, whereas the dashed red one 
includes also contributions from E + and E° production, as- 
suming cr(E ) = <t(E + ) = 0.021 p,b, <ti.826(A) = 7.9 and 
o"i.775(A) = 6.1 ph. The upper magenta dotted histogram 
represents a simulation where the E + total cross section is 
taken to be six times that of E°, which is the limit allowed 
by the triangle inequality at e = 129 MeV. 
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FIG. 5: (Color online) Laboratory differential cross sections 
for inclusive K + production at the five measured energies, 
evaluated over the angular region -&k < 4°. The overall un- 
certainty of 7% coming from the luminosity has not been in- 
cluded. The experimental results are compared with Monte 
Carlo simulations; A production (black dotted line), A + E° 
production (blue dashed line), A + E° + E + production (red 
solid line). 



The average momentum px of the K + meson, detected 
in each of the telescopes, as well as the corresponding 
values of the inclusive double-differential cross sections, 



evaluated using Eq. ©, are presented in Table UTTl These 
have then been summed and averaged over the kaon an- 



8 



gular range $#-+<4° to give the inclusive cross sections 
shown in Fig. [5] for the five beam energies studied. The 
data are compared with simulations of A, E° and S + pro- 
duction discussed in Sec. [TTl where the individual contri- 
butions from the three reaction channels have been nor- 
malized to the total cross sections deduced from the K + p 
and K + ir + correlation data, as described in Sec. [V] Al- 
though there is fair agreement between the simulations 
and the data for all five energies, it has to be stressed 
there is considerable ambiguity in the dominant A pro- 
duction cross section because in this case the reaction is 
measured far above threshold. The COSY-TOF results 
show that there are then strong deviations from phase 
space for both the angle and momentum of the K + 0. 



C. Kaon-proton coincidence measurements 

The K + p missing-mass spectra measured at four beam 
energies are presented in Fig. |6] Due to technical prob- 
lems in the preparation of the K + p correlation trigger, 
the spectrum measured at 1.826 GeV could not be recon- 
structed reliably and is therefore not shown. However, 
this problem does not affect either the K + inclusive or 
the K + 7T + correlation data collected at this energy. 

Figure [6] shows clear peaks associated with the pp —> 
K + pA and pp — > K + pYP production processes, where 
the "p" corresponds to the proton that is detected in 
ANKE. For events that lie within the E° peak, the distri- 
butions of K + missing masses are consistent with predic- 
tions based upon a phase-space production model, with 
no evidence for any strong YPp final state interaction. 

In addition to the direct protons, there are large "in- 
direct" physics backgrounds, where the K + is measured 
in coincidence with a proton that arises from the decay 
of one of the hyperons produced. The values of the rele- 
vant decay fractions are given in Table H Both types of 
contribution from the three reaction channels have been 
simulated and traced through the GEANT4 model of the 
ANKE spectrometer. In view of their proximity to the 
thresholds, as well as the evidence discussed above, the 
contributions to the spectra from the two £ channels 
have been estimated using simple phase-space descrip- 
tions. However, A hyperon production, which dominates 
the distributions, has been simulated using the different 
models mentioned in Sec.|TTJ In all cases the protons that 
were produced in the decay of the three hyperons were 
assumed to be distributed according to phase space. 

The numbers of detected A and E° events could be 
fixed reliably from fits to the experimental peaks and 
do not depend significantly on the model used for the A 
production. The indirect protons from A decay lead to 
a missing-mass distribution that dies off well before the 
kinematic limit, as can be judged from the simulations 
shown in Fig. [6j Above the A limit there are contribu- 
tions coming from protons from E° and E + decay and 
these give rise to distributions with somewhat different 
shapes. The numbers of E° events expected in this region 




mm(Kp) [GeV/c 2 ] 

FIG. 6: (Color online) The K + p missing-mass spectra mea- 
sured at four beam energies (closed symbols) compared to 
simulations normalized to the total cross sections listed in 
Tables [V] and IVII In addition to the protons from pp — ► 
K + pA/T,° , these include contributions from indirect K + p cor- 
relations involving decay protons arising from the production 
of the A (histogram 1), E° (histogram 2 (blue)), and E + (his- 
togram 3 (red)). As explained in the text, the K + p correla- 
tion data were unreliable at 1.826 GeV and so results are not 
shown for this energy. 



can be estimated from the population in the TP peak and, 
by subtraction, the numbers of E + events can be esti- 
mated. The resulting values for the E + production cross 
section are 0.25 ± 0.05, 0.41 ± 0.09, and 1.09 ± 0.16 ^b at 
1.920, 1.958 and 2.020 GeV, respectively. Only statisti- 
cal error bars are quoted and it is hard to obtain robust 
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Pk 
MeV/c 




da / dQdp 
nb/(srMeV/ 


=) 




#(1.826/1.775) 


231 


4.7 ±0.5 


4.3 ±0.5 


4.2 ±0.5 


6.1 ±0.7 


7.5 ±0.9 


0.91 ±0.02 


258 


6.1 ±0.7 


6.1 ±0.7 


8.1 ±0.9 


11.8 ± 1.3 


12.0 ± 1.3 


0.98 ±0.02 


285 


7.6 ±0.9 


7.1 ±0.8 


11.0 ± 1.2 


12.8 ± 1.4 


12.6 ± 1.4 


0.94 ± 0.02 


315 


8.6 ± 1.4 


9.4 ± 1.3 


14.9 ± 2.1 


15.8 ± 2.4 


16.7 ±2.5 


0.90 ±0.04 


348 


12.6 ± 1.4 


13.6 ± 1.4 


17.6 ± 1.9 


19.2 ± 2.1 


18.9 ± 2.0 


1.03 ±0.02 


381 


14.4 ± 1.4 


14.2 ± 1.3 


19.7 ± 1.8 


21.3 ± 2.0 


19.4 ± 1.9 


0.93 ± 0.02 


415 


17.1 ± 1.6 


16.4 ± 1.5 


22.4 ± 2.0 


22.4 ± 2.1 


25.0 ± 2.3 


0.95 ± 0.02 


452 


18.9 ± 1.7 


18.7 ± 1.7 


22.3 ± 2.0 


26.7 ±2.5 


27.1 ± 2.5 


1.00 ±0.02 


488 


24.6 ± 2.2 


23.6 ± 2.1 


28.3 ± 2.5 


30.1 ± 2.7 


34.0 ±3.0 


0.99 ±0.02 


525 


27.2 ± 3.3 


28.4 ±3.3 


30.7 ±3.7 


32.6 ±3.9 


38.3 ±4.6 


1.01 ± 0.02 


559 


27.8 ±3.6 


29.7 ±3.7 


34.7 ±4.4 


37.4 ±4.8 


41.4 ±5.3 


1.01 ± 0.02 


595 


29.0 ±4.0 


30.6 ±4.1 


34.2 ±4.6 


36.6 ±5.1 


43.2 ± 5.9 


0.99 ± 0.02 


T p [GeV] 


1.775 


1.826 


1.920 


1.958 


2.020 





TABLE III: The K + laboratory double-differential cross sections measured at five beam energies are integrated over $ K + < 4° 
and momentum bins of width ±12 MeV/c. The errors do not include the systematic uncertainty of 7% coming from the 
normalization. The final column shows the ratio of K + count rates measured at 1.826 and 1.775 GeV in the interval i9 K + < 12° 
and ±12 MeV/c momentum bins. 



estimates of the systematic uncertainties due to the prox- 
imity of the kinematic limit. Conservative upper limits 
on the S + production cross section can be derived if one 
assumes that the whole contribution to the spectra above 
the A limit comes from S + production. This extreme 
hypothesis gives values of 0.51, 0.83, and 1.62 pb at the 
three energies. These results are given in Table IVT1 along 
with those obtained from the K + ir + analysis. 

The histograms in Fig. [5] present our estimations for 
all three reaction channels, normalized to the total cross 
sections listed in Tables [V] and IVII If the value of the 
E + cross section at 1.958 GeV were as large as that ad- 
vocated by the COSY-11 collaboration [8j, this would 
require boosting the red histogram in Fig. [6] by about 
two orders of magnitude. 



D. Kaon-pion coincidence measurements 

Below the threshold for the pp — > K + nAir + reaction 
(Tthr = 1.975 GeV), the K+tt + coincidences arising from 
the hyperon decay S + — > mr + are a unique signature 
of the pp — > K + nTi + reaction. The background from the 
K + nAn + channel will remain negligible also at 2.02 GeV 
because of the very small phase space and the fact that at 
the much higher energy of 2.88 GeV it comprises only 4% 
of the E + production cross section @. The smallness of 
the K + nAn + production rate was also confirmed by our 
studies at 2.16 GeV, where no evidence for this channel 
was found in the kaon-proton and kaon-pion coincidence 
spectra 0]. 

The measured kaon and pion momentum spectra from 
K + n + coincidence events are presented in Fig. where 
they are compared with phase-space simulations of the 
pp — > K + nH + reaction. The number of random coin- 
cidences that survive the analysis procedure discussed 
in Sec. IIIII can be judged from the spectra measured at 



T p = 1.775 GeV, which is below the £+ threshold. The 
kinematically allowed region in the two dimensional plane 
of the K + and 7r + momenta has been determined from 
simulations carried out at each energy. The resulting 
kinematical cuts have been applied, not only on the dis- 
tributions above threshold, but also on the 1.775 GeV 
data. To compensate for the random background, the 
number of events estimated from the subthreshold data 
has been subtracted from the data at the different ener- 
gies, taking into account the relative total luminosities. 



e 


Ltat 


N tot 


acc 


MeV 


Pb" 1 


event 


% 


13 


3.53 


16 


2.3 


47 


2.17 


54 


0.6 


60 


1.49 


54 


0.5 


82 


1.40 


69 


0.1 



TABLE IV: Input information for the E + total cross sec- 
tion estimation on the basis of the K + ir + correlation data. 
Values are given for the excess energy (e), total luminosity 
(Ltot), accumulated number of K + n + coincidences after back- 
ground subtraction (N to t), and total geometrical acceptance 
used (e^+). 

After subtracting the background bin-by-bin and cor- 
recting for efficiencies, the K + and tt + momentum spec- 
tra measured above the S threshold are compared in 
Fig. [7] with phase-space simulations. Although the num- 
bers of events are small and the fluctuations large, both 
the K + and 7r + distributions are in reasonable agree- 
ment with the simulations, within the uncertainties of 
the background subtraction procedure. The number of 
K + 7r + coincidences extracted from the spectra, together 
with estimates of the overall geometrical acceptance (see 
Tab. lIVp . including all the efficiencies, lead to evalua- 
tions of the total cross sections on the basis of Eq. ([5]). 
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FIG. 7: Spectra of the kaon and pion momenta from identified 
K + tv + pairs measured at four energies above the E + threshold 
and one below. The pions are taken over the full horizontal ac- 
ceptance of ANKE (« 15°) whereas the kaons are cut at 12°. 
Random background has been subtracted bin-by-bin from the 
above-threshold data by scaling the 1.775 GeV results by the 
relative luminosity. The spectra are compared to phase-space 
simulations of the pp — ► K + nY, + reaction channel. 



These are presented in Table IVI1 where the very satisfac- 
tory agreement with the results derived from the K + p 
correlation data is demonstrated. 



V. TOTAL CROSS SECTIONS FOR HYPERON 
PRODUCTION 

The results obtained for the total cross sections of 
hyperon production in proton-proton collisions are pre- 
sented in Table [V] and ED The values in the A and E° 
cases are extracted from the fits to the two peaks in the 



K + p correlation spectra of Fig. [6] Since the data were 
taken far from the A threshold where the ANKE accep- 
tance is low, the uncertainties here are dominated by the 
ambiguities in the models used to describe the differential 
distributions. Any model dependence is far weaker for E° 
production and the contributions from other systematic 
effects, including luminosity, are of a similar size. 
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1.775 


70 ±1 


6.1 


2 


18 


A 


1.920 


122 ± 1 


16.9 


2 


21 




1.958 


135 ± 1 


16.1 


2 


21 




2.020 


157 ±1 


18.0 


2 


31 




1.920 


45 ± 1 


0.24 


5 


14 


E° 


1.958 


58 ±1 


0.51 


5 


14 




2.020 


80 ±1 


1.30 


5 


14 



TABLE V: Total cross sections for A and S° production mea- 
sured in proton-proton collisions. 
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0~Kp 
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GeV 


MeV 
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1.826 


13 ± 1 


0.011 


38 


17 






0.045 


1.920 


47 ±1 


0.20 


16 


17 


0.25 


21 


0.51 


1.958 


60 ±1 


0.39 


15 


17 


0.41 


23 


0.83 


2.020 


82 ± 1 


0.76 


14 


17 


1.09 


15 


1.62 



TABLE VI: Total cross sections for pp — » A" + nE + produc- 
tion measured in this experiment. The values in the column 
marked axn were obtained from the study of the K + n + coin- 
cidence data whereas those in column <jk p followed from the 
study of the high K + p missing-mass spectra of Fig. [6] As 
discussed in the text, it is very hard to quantify satisfacto- 
rily the systematic uncertainties in this case. The upper limit 
quoted in the last column at 1.826 GeV was derived from the 
study of inclusive K + production above and below threshold 
shown in Fig. [4] At the higher energies, the upper limits were 
obtained by assuming that all the events in Fig. [6] that fell 
above the kinematic limit for the A were associated with S + 
production. 

Table fVTl shows the values of the S + production cross 
section extracted from both the K + ir + and K + p cor- 
relation data shown in Figs. [7] and |H1 respectively. As 
previously remarked, experimental difficulties at T p — 
1.826 GeV made the K + p correlation data unreliable at 
this energy, though this had no influence on the K + ir + 
measurement. The principal errors in the K + n + case 
arise mainly from the low statistics. These are much 
better in the K + p coincidence measurement but the dif- 
ficulties associated with evaluating the acceptance at the 
kinematic limit makes it hard to quantify the system- 
atic errors. Despite these uncertainties, the agreement 
between the two sets of measurements is very reassur- 
ing. Upper limits on the E + production cross section are 
also quoted in Table IVII That at the lowest energy was 
extracted at the 98% confidence level from the compar- 
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ison of the inclusive K + production data at 1.826 and 
1.775 GeV shown in Fig. 2) At the three higher energies 
the upper bounds were determined from the assumption 
that all the events above the A limit in the K + p missing- 
mass spectra of Fig. [6] corresponded to E + production. 
At 1.958 GeV this is about a factor of 50 below the re- 
sult of the COSY-11 collaboration Q. It can be seen from 
the tables that the E + production cross sections are only 
a little less than those of E° but that the difference is 
within the error bars. 

The results for all three hyperons are presented in 
graphical form in Fig. [5] It is seen from this that our 
values of the production cross sections for both A and S° 
are consistent with the world data. However, it has to 
be stressed once again that a large systematic error bar 
had to be ascribed to the A value because of the small 
ANKE acceptance. As is well known, it is important to 
introduce the Ap final state interaction to describe the 
energy dependence of the pp — > K + pA total cross sec- 
tion. On the other hand it is seen that the pp — > K + pTP 
variation is much closer to phase space, though our data 
might give a possible hint of some small repulsive FSI. 

The E + production cross section in pp collisions shows 
no sign of any threshold enhancement of a type that could 
explain the large values found by the COSY-11 collabo- 
ration Q and, indeed, their values are already completely 
excluded by the upper bounds provided by the K + inclu- 
sive and K + p correlation data of Figs. 21 [5] and [5] The 
theoretical models [HI, [33| whose predictions are shown 
in Fig. [SJd were motivated by the COSY-11 data and fail 
completely to reproduce our results. 

The ratio of S + to S° production seems to lie between 
about a half and one, though the difference is possibly 
within the error bars. The energy dependence shown 
in Fig. [8)3 seems a little steeper than that suggested on 
phase-space grounds and, if this is so, it could be an indi- 
cation for a damping effect near threshold due to a S + n 
final state interaction. There is ample evidence from the 
study of the threshold cusp in K~d — > n~Ap 34| that the 
Ap and £A channels are strongly coupled. The enhance- 
ment in the K + missing mass in the inclusive pp — > K + X 
data at the £ N threshold [ill, Ell is a possible sign that 
flux is being removed from the pp — » K + NT, channels in 
this region. 



VI. CONCLUSIONS 

We have reported measurements of the pp —> K + nT, + 
reaction at four energies close to threshold using two 
techniques. Instead of identifying the neutron, as was 
attempted at COSY-11 Q, the K + was detected in coin- 
cidence with the 7r + that came from the decay of the E + 
hyperon. Below the threshold for K + nAir + production, 
this decay is the only source of such correlations. The 
numbers of random K + n + coincidences were small and 
could be estimated from data taken below the threshold 
for E + production. These results were confirmed within 
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FIG. 8: (Color online) Total cross sections for hyperon pro- 
duction in proton-proton collisions as functions of excess en- 
ergy obtained in this experiment and compared to published 
data. Panel (a): Results for A (black stars) and E° (blue 
stars). The data points obtained in other experiments at 
COSY and other facilities are shown by open and closed sym- 
bols (circles for A and squares for E°) [tHl Hfl, 1, Gl, ETj, . 
The lines represent parameterizations of the world data from 
Ref. [14| . Panel (b): E + production. Closed stars represent 
the current data set obtained from the K + ir + coincidence 
spectra. The open star corresponds to our previous measure- 
ment [3]. The conservative upper limits given in Table IVTI 
are shown by the four short horizontal lines. The COSY- 
11 results [S| are shown by open triangles while the closed 
ones are taken from a compilation of higher energy data [3ll ] . 
The dotted line corresponds to a pure phase-space variation 
with a normalization chosen to pass through our central en- 
ergies. The solid line is a parameterization that includes a 
strong nE + FSI and energy dependent matrix element [3^ ]. 
The broken line is taken from the microscopic calculation of 
Ref. M. 



error bars from the study of the high end of the K + p 
missing-mass spectra, though it has to be stressed that 
systematic effects are harder to control in this measure- 
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ment. Nevertheless, conservative upper limits on E + pro- 
duction could be obtained by assuming that the S° events 
did not populate the high K + p missing-mass distribu- 
tions. In complete contrast to the COSY-11 results Q, 
we found cross sections for the pp — » K + nS + reaction to 
be even a little smaller than those of pp — > K + pYP . 

The large COSY-11 E + production rates could also be 
excluded by the study of inclusive K + production [l2j]. 
However, this does not provide the best estimates of the 
pp — » K + nYi + cross section because the K + spectrum is 
dominated by the K + pA channel, unless S + production 
is anomalously large near threshold §. The uncertainty 
in the £+ cross section deduced from such data cannot 
therefore be smaller than the sum of the uncertainties in 
the models describing A and E° production, which is why 
these data only provide upper limits on E + production. 

The data on A production were taken far above thresh- 
old where the acceptance of ANKE for the pp — > K + pA 
reaction is small. The predicted numbers of inclusive K + 
depend sensitively upon the variation of the production 
with angle and momentum. In addition to the strong 
distortion due the Ap final state attraction, the COSY- 
TOF data @ show significant angular variation. As a 
consequence, one cannot extract total cross sections for 
A production in a model-independent way from an ap- 
paratus with a small acceptance in angle. The situation 
would be even more uncertain if the strong Ap ^ EiV 
channel coupling [34j gave rise to a rapid K + momentum 
dependence. Our aim was merely to show that, with rea- 
sonable model assumptions, we could achieve total cross 
sections that were not inconsistent with world data. 

The experimental situation regarding acceptance is far 
more comfortable for the K + pT,° channel in this respect 
because the results were obtained not far above thresh- 
old. Furthermore, it is believed that the YPp final state 
interaction is much weaker than that of the Ap and 
our data agree with this. The resulting values of the 
pp — > K + pTP cross section shown in Fig. (5£i are com- 
pletely consistent within error bars with published data 
and this supports the validity of our experiment and its 
analysis. 

Although the K + tt + coincidence method leads to only 
small numbers of events within our acceptance, these are 
very clean and so the total cross sections presented in 
Fig. [Bt> have fewer systematic uncertainties than those 
coming from the K + p correlation data. It may there- 
fore be interesting to note that, taken together with our 
2.16 GeV point [7], the data seem to indicate a some- 
what steeper rise than expected on the basis of phase 
space. This is definitely the case if the data have to join 
smoothly on to the bubble chamber measurements shown 
in Fig. E 



Attractive final state interactions, such as those in the 
pp or Ap systems, give rise to a very strong energy de- 
pendence because they can generate virtual state poles 
in the production amplitudes. This is not the case when 
there is repulsion between two final particles and so this 
would, in general, lead to much smaller effects. There is, 
however, a third type of final state interaction, namely an 
absorptive one where the final wave function is damped 
at small distances rather than being repelled. In the 
present case this could arise from the coupling of the UN 
and the Ap channels, which is known to be very strong at 
low energies [13]. It is therefore possible that flux from 
the K + NYi channels is lost to feed the K + pA channel. 

If the coupled channel final state interaction is signifi- 
cant, it would only be relevant in the isospin I = \ state 
and this would have twice as big an effect on K + nT, + 
production compared to K + pYP. It is perhaps note- 
worthy that the parameterization of the world data for 
pp — > K + pYP shown in Fig. [8^, do not seem to indicate 
as steep a rise as for pp — » if + n£ + , though this is not 
the case for our own S° data. 

As stressed in Sec. HI1 there are isospin links between 
the amplitudes for pp — > K + pYP, pp — > K + nY, + 1 and 
pp — » K Q pY. + , but these all refer to an initial N N isospin- 
one state. The isospin-zero amplitude for pn — » K + pY<~ 
is independent of these and so extra help for the modeling 
might be provided by measurements in this channel. This 
will be possible at COSY-ANKE through the study of 
quasi-free production on the deuteron with the detection 
of the spectator proton 36]. 

The existing theoretical models for pp — > K + nT, + 
production [32|, bear no relation to our results. If 
the reaction were driven by the exchange of non-strange 
mesons exciting a A isobar (33| , this would strongly favor 
the production of the S + compared to the S°, which is 
in conflict with the data. It would also not contribute 
to the isospin-zero production of Further theoreti- 
cal work is clearly needed to understand the full data set. 
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